SEARCH FOR WMAP- COMPATIBLE SIMPLE 50(10) SUSY GUTs 



SEZEN SEKMEN 

Department of Physics, Middle East Technical University, 
TR- 06531, Ankara, Turkey 

Unification of GUT-scale t — b — r Yukawa couplings is a significant feature of simple SO{10) 
SUSY GUTs. Here we present the results of a search that used the Markov Chain Monte 
Carlo technique to investigate regions of Yukawa unification and WMAP-compatible dark 
matter relic density in 5*0 (lO)-like MSSM parameter spaces. We mention the possible LHC 
signatures of Yukawa unified scenarios and discuss the consequences for dark matter. 

1 Introduction 

Grand Unification is regarded as an inspirational ingredient of models that claim to explain the 
fundamental laws of nature. A highly motivated scenario in this context originates from grand 
unification via the SO{10) gauge group ^ . Simple supersymmetric implementations of SO{10) 
GUTs unify all matter fields in each generation within a 16-dimensional irreducible representa- 
tion and two Higgs doublets of the MSSM within a 10-dimensional irreducible representation. 
Such a formalism automatically includes heavy right-handed neutrino states and the resulting 
structure of the neutrino sector implies a successful theory of baryogenesis via intermediate scale 
leptogenesis. Moreover the SO{10) models are left-right symmetric and this enables them to 
provide a solution to the strong CP problem and to naturally induce R-parity conservation. 

Besides gauge coupling unification, 50(10) SUSY GUTs additionally require the unification 
of 3rd generation Yukawa couplings at the GUT scale {Mqut)- This is explicitly seen from the 
expression of the superpotential above Mqut^ which takes the form / 3 /'^le'^ie^io + • • •• An 
exact unification occurs at tree level while several percent corrections arise at the loop level. As 
a result we can assume that any sign of Yukawa unification from observations could be a hint 
to the existence of 50(10) SUSY GUTs. 

Our aim in the study is to investigate the characteristics that arise in a SUSY model when 
GUT scale Yukawa unification is imposed and to determine the experimental signatures that 
would distinguish such Yukawa- unified models from the others. In this context we assume a 
theoretical framework where nature is explained by an 50(10) symmetry above Mqut- Then 
at MguTi 50(10) breaks to MSSM plus some heavy right-handed neutrino states. At the weak 
scale, the content of the theory is equivalent to that of the MSSM. 

The GUT scale soft SUSY breaking parameters are constrained by the requirement of the 
50(10) symmetry. Unified representations would favor common SSB masses "mie" for the 
scalars and "mio" for the Higgses, but in order to achieve REWSB, SSB Higgs masses should be 
split, satisyfing ttih^ > Here we examine two different methods to generate the necessary 

Higgs splitting: The first approach defines the Higgs masses as ^'h^ ^ = T 2M|). Here 
splitting is parametrized by Mjj, which is the magnitude of the D-terms in the scalar potential 



of the extra U{1] group that is a by-product of the 50(10) breaking. The parameters of this 
GUT scale Higgs input (GSH) scenario are 



mi6, mio, M|), mi/2, Aq, tan/3, sgn{ii) (1) 

The socond approach was put forward in order to generate Yukawa unified solutions with low fi 
parameter and low m^. Such solutions were found to exist by Blaszek, Dermisek and Raby at 
a study where they assumed perfect Yukawa unification at Mqut and made a fit to the weak 
scale observables ^ . In order to seek similar solutions, we start with GSH parameters at GUT 
scale, but additionally provide ji and tua as inputs. We run tuh^ ^ down, and si Q — MsuSY 
we compute what m/f^, tuh^ should have been in order to give our input ji and m^, and run 
back up using these new boundary conditions. This weak scale Higgs (WSH) scenario has the 
parameters 

mi6, mio, M|), mi/2, tan/3, m^, (2) 

We take the GSH and WSH scenarios and search in their parameter spaces for regions 
having a good Yukawa unification where Yukawa unification is parametrized hy R— • 
Additionally we seek sub-regions that are consistent with the WMAP measurements of dark 
matter relic density Q.h?^. There have been previous searches using the GSH input based on 
random scans and they were able to achieve less than few percent of Yukawa unification for 
> 0^. However the dark matter relic densities for these solutions were always much higher 
than the WAMP upper bound. Here we implement the Markov Chain Monte Carlo (MCMC) 
technique which enables a much more efficient scanning of multi-dimensional parameter spaces. 
The following two sections summarize the MCMC technique, the characteristics of the regions 
found by utilizing it, and dark matter-related consequences of the 50(10) SUSY GUTs. 



2 The MCMC Search and the Yukawa- unified Solutions 

A Markov Chain is a discrete time, random process where given the present state, the future 
state only depends on the present state, but not on the past states^. The MCMC samples from 
a given parameter space as follows: It takes a starting point, and it generates a candidate point 
from the starting point using a proposal density Q{x^] x^). The candidate point is accepted 
to be the next state x'^^^ if the ratio p — p(^tjg(^c!^tj (where P{x) is the probability calculated 
for the point x) is greater than a uniform random number a — U(0, 1). If the candidate is not 
accepted, the present point x^ is retained and a new candidate point is generated. By repeating 
this procedure continuously the Markov Chain eventually converges at a target distribution 
around a point with the highest probability. 

Our MCMCs were directed to approach regions with R - 1.0 and 0.094 < Qh^ < 0.136. Here 
we used a Gaussian distribution for the proposal density Q, and approximated the likelihood of 

2 

a state to e~^^'^. We also chose multiple starting points (~ 10) in order to search a wider range 
of the parameter space. We used ISAJET 7.75 for sparticle mass computations and micrOMEGAs 
2.0.7 for DM relic density calculations. The MCMCs successfully located some regions with 
good R and fi/i^. Here we will mostly emphasize the results from the GSH scenario. 

Figure [T] shows the compatible regions in tuiq — mio and uiiq — mie/^o planes. The light 
blue dots have R < 1.1, the dark blue dots have R < 1.05, the orange dots have i? < 1.1 
plus Qh^ < 0.136 and the red dots have i? < 1.05 plus Qh^ < 0.136. We see that Yukawa 
unification occurs only at the regions where the input parameters are strongly correlated, having 
mio — 1.2mi6 and Aq 2^ —(2 — 2.1)mi6. A good DM relic density is achieved only at the 
constrained regions that have mie ^3 — 4 TeV. Further search showed that mi/2 takes the 
lowest possible values for a given mie, generally giving ~ 100 GeV, and decreases steadily with 
increasing tuiq. 




Figure 1: Yukawa- unified GSH points found by MCMC in the mio vs. mie plane (left) and the Ao/mio vs. 
mi6 plane (right); the light-blue (dark-blue) points have R < 1.1 (1.05), while the orange (red) points have 

R < 1.1 (1.05) plus Qh^ < 0.136. 
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Figure 2: Yukawa- unified GSH points found by MCMC in the vs rrig plane (left) and mi^o vs. m^o — m^o 

plane (right); color code as in Figure [T] 



These highly confined parameter regions lead to strongly constrained mass spectra, and 
hence to significant LHC signatures. We see that Yukawa-unified solutions are distinguished 
by their heavy lst/2nd generation scalars (> 2 TeV), lighter 3rd generation scalars (~TeV) 
and light gauginos (few hundred GeV). All Higgses except are about 1 — 3 TeV. Figure [2] 
shows the distribution of selected points on m^^ vs rrig plane (left), and on m^o vs m~o — m~o 
plane (right) for the GSH scenario. The requirement of Qh^ < 0.136 favors a gluino mass range 
around 350 — 450 GeV, which means we would expect a large amount of gluino pair production 
at the LHC with cross sections about ^ 100 pb. The gluinos decay via 3-body channels such as 
g Xi^^i X2^^5 X^tb/bt^ since 2-body channels are closed due to the high squark masses. On the 
other hand favored X2 — xt mass range is 100 — 150 GeV, which leads to gaugino pair production 
cross sections about ^ 10 pb, while m~o ^ 50 — 75 GeV. The preferred mass difference m~o — m~o 
is 52 — 65 GeV which is smaller than ruz^h^^ therefore X2 decays are dominated again by 3-body 
channels as X2 ^ ^^Xi, QQXi, l^Xv 

As a result we expect the 50(10) models to manifest themselves as multi b-jet plus low 
missing Et final states at the LHC. Additionally it would be possible to investigate the OS/SF 
dilepton channels where the dilepton invariant mass is bound by m ~o —^^0 and try to reconstruct 
the g bbx2 ~^ bbllxi cascades. 
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Figure 3: Yukawa- unified points found by MCMC on the ttia — 2m~o vs. m^o — 2m~o plane for GSH (left) nd 

WSH (right) scenarios; color code as in Figure [1] 

3 Consequences for Dark Matter 

The majority of solutions shown in Figures [1] and [2] have excess DM rehc density, while only 
a small portion of them gives Qh'^ < 0.136. To investigate the mechanism that provides the 
efficient annihilation in the DM-allowed solutions, we check the behavior of solutions in m^o 
and rriA mass resonances. Figure [3] shows the distribution of Yukawa-unified points on the 
rriA — 'im^o vs. m^o — 2m~o plane for GSH (left) and WSH (right) scenarios. In the GSH plot all 
DM-allowed solutions are on the m/^o 2^ 2m^o line, which shows that the relic density is reduced 
by annihilation via a light Higgs resonance. On the other hand tua > 2m ~o, so there are no A 
resonance solutions. Turning to the WSH scenario we see that annihilation via both and A 
resonances are at work. Actually the majority of solutions are generated by the latter due to 
the relatively small A masses allowed within the WSH scenario. However all of these solutions 
have Bg jJ^jJ^ branching ratios higher than the latest reported CDF upper limit 5.8 x 10~^. So 
these A resonance points are ruled out, leaving us with only the resonance solutions. 

One could also devise alternative methods for reducing the excess DM relic density. One way 
could be to assume that Xi is not the LSP, but can decay to other candidates such as gravitino 
or axino via the mode Xi ~^ jG/a. Lifetime of Xi would be long enough to let it escape the 
detectors. The resulting relic density would be ^ = (m^ ~/m~o)fi^o since the Gs/as inherit 

the thermally produced neutralino relic number density. G LSP can only reduce the relic density 
a few times, which is not satisfactory for our case, but axinos with < 1 MeV would allow for 
a mixed cold/warm DM solution which can reduce the relic density below the WMAP bound ^. 

Another method to reconcile ^h? is to relax some universalities in the GUT scale SSB terms. 
For example, increasing the U{1) gaugino mass term Mi (while keeping M2^3 = '^1/2) brings m~o 
close to m~± , hence making x? more wino-like and inducing bino-wino coannihilation. A further 
possibility is to lower the lst/2nd generation masses mi6(16) (while keeping mi6(3) = ^le)? 
which enables neutralinos to annihilate via light qR exchange and leads to neutralino-squark 
coannihilation. 

4 Conclusions 

By performing scans on the parameter space of simple 50(10) SUSY GUT scenarios using the 
MCMC technique, we showed that solutions with both 5 — 10% Yukawa unification and WMAP- 
compatible ^h? can exist around mie ~ 3 — 4 TeV. These regions defined by strictly constrained 
relations among the GUT scale inputs generate special sparticle mass relations that lead to 
distinguishable signatures at the LHC. With multi-TeV scalars, 350 — 450 GeV gluinos and 



50 — 150 GeV light gauginos, we expect dominant gluino production followed by 3-body cascade 
decays which will end up in b-rich multijet final states, occasionally including OS/SF lepton 
pairs from X2 decays. Moreover the possibility to lower ^Ih^ by assuming a LSP or introducing 
SSB non-universalities marks wider parameter space regions as compatible. So we can conclude 
that 50(10) SUSY GUTs provide motivated scenarios with robust signatures relevant to be 
tested soon at the turn-on of the LHC. 
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